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Abstract 

Quantum computation is based on implementing selected unitary transformations which 
represent algorithms. A generalized optimal control theory is used to find the driving field that 
generates a prespecified unitary transformation. The approach is illustrated in the implementation 
of one and two qubits gates in model molecular systems. 
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A universal model of a quantum computer can be constructed from an array of two level 
systems used as registers (qubits). Any general quantum gate can be decomposed to a one 
and two qubit unitary transformation [[!]]. Thus a physical realization of a quantum computer 
should be able control such unitary transformations by use of external driving fields. Real- 
izations based on nuclear magnetic resonance techniques already have been demonstrated for 
non-trivial quantum gates. For example, a three qubit quantum Fourier transform has been 
performed by a sequence of pulses corresponding to one and two qubit transformations [§. 
This implementation was based on addressing each qubit independently by their spectral 
separation. However, when many qubits become involved in a computation, the spectrum 
becomes congested so that it becomes more difficult to address each qubit individually. The 
total fidelity of the algorithm will depend on the accumulation of errors at each step. An- 
other difficulty results from decoherence processes which are unavoidable due to coupling to 
the external environment. These processes will degrade the performance in proportion to 
the time required to carry out the computation task. Thus it is desirable to minimize the 
number of computation steps and the total computation time. 

An alternative possibility for implementing quantum computing is the use of molecules 
driven by shaped light pulses. A hypothetical model can be based on using the vibrational 
and rotational states as registers. A shaped light pulse can write data as amplitudes on these 
states. The computing algorithm is then implemented by performing a unitary transforma- 
tion employing a second shaped pulse. In the end the output can be read by a probe pulse. 
This approach implies changing the computation model by tying together many single and 
two qubit operations into one resulting in a combined unitary transformation. 

Such an approach has been used for the experimental implementation of elements of co- 
herent computation in Li2 H f|. For example, two states of the ro- vibrational manifold 
of an electronic state can form a single qubit. By taking advantage of multiple electronic 
states, the evolution of the molecular system can be controlled using electromagnetic fields 
in the optical region realized by shaped laser pulses. The fast development of pulse shaping 
technology can make possible the implementation of simple quantum gates in molecules. 
This program has to overcome interference from the large number of other molecular levels 
coupled to the field but not assigned to the qubits. The task becomes therefore to imple- 
ment the quantum algorithm on the molecular levels used as registers while avoiding the 
intervention of other states from the same system. 
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Tesh et. al. f5|, proposed the use of optimal control theory (OCT) to calculate the 
field which can induce a specific transformation used for quantum computation. Originally, 
OCT was designed as a method to obtain a light field which could induce a specific state 
to state transformation [7|. In the quantum computing context, the goal is to obtain the 
optimal pulse that induces a given unitary transformation, irrespective of the initial state 
of the system. The present approach generalized OCT to obtain directly the driving field 
that induces a target unitary transformation in the system. The approach is based on the 
equation of motion of the unitary transformation by making use of the main ingredients of 
OCT. This allows the utilization of the large number of tools that have been developed in 
this field. 

The model system consists of a free Hamiltonian H controlled by an external field e(i), 

H(t) = H -Mt). (1) 

where p, is a system operator. In the molecular system, fi is the transition dipole operator 
and e(t) describes a shaped short light pulse. For simplicity, the field e(t) is assumed real but 
the generalization to an electromagnetic field with two independently controlled polarizations 
|| is straightforward. 

The algorithm is represented in the target time T by a unitary transformation U(T) 
generated from the Hamiltonian Eq. (|]); 

^ --><*). ( 2 , 

with the initial condition TJ(t = 0) = 1, where 1 denotes the identity operator. The objective 
is to obtain the optimal driving field e(t) that induces a given unitary transformation O at 
t = T, i.e., U(T) = e^O. (ft denotes a physically irrelevant global phase that is related to 
the energy origin and is therefore uncontrollable by the field. 

The task is a typical inversion problem which can be solved by employing a variational 
procedure maximizing the projection of the generated operator on the target operator, 

|r| = \Tr{6 j V(T)}\ , (3) 

where the projection (A • B) is defined by Tr{A^B}. The functional r is a complex number 
inside the circle \r\ < Njj, where Nh is the dimension of the Hilbert space of the system. 
Equality is reached only when the argument of the trace is e 1 ^ 1 and then a maximum 



of |r| is equivalent to U(T) = e^O. The optimal solution is then found by maximizing 
the functional Eq. (D), with respect to the control field e(t). Since a direct algorithm to 
maximize |r| was not found, a working alternative is used based on formulating the problem 
as the optimization of Re[r], or of Im[r], or of a linear combination of both. For simplicity, 
the optimization of the real part represented by the functional J = Re[Tr{6 U(T)}] is 
considered. Two constraints are introduced || [7[], the first restricts the dynamics to obey 
the Schrodinger equation Eq. (|2]), and the second restricts the total field energy. Using 
Lagrange multipliers a modified functional is obtained, 



J = Re 
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A / -TT-Mt^dt , (4) 
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where B(t) is an operator Lagrange multiplier, A grange multiplier and s(t) is 

a shape function which turns the pulse on and off || . The use of more elaborate constraints 
and choices of A, allows a higher degree of control on the shape of the optimal pulse [[IIJ [TT 



Applying the calculus of variations, 5 J = with respect to B, U, and e, a set of equations 
is obtained: a) The Schrodinger equation Eq. (0) with the initial condition XJ(t = 0) = 1 
for U; b) The inverse Schrodinger equation 

^4b W H W (5) 

with the condition B(t = T) = for B; c) The field equation: 

e (t) = -|^Im[Tr{B(t)AU(t)}], (6) 

Eq. (0) and (||) represent two counter currents with information from the initial condition 
and the target unitary transformation respectively. The equations are solved iteratively, the 



Krotov method p2| , similar to the methods described in Ref. ||13|| , was found to be the 
most efficient. The input is a "guess" field, e^°'(t), so that in the k iteration [k = 1,2...): 
(i) B*' \t) is propagated backwards from t = T to t = using Eq. (||) and e^" 1 -*; and (ii) 
U (t) is propagated forward using Eq. (||) and e^ k > is evaluated using 

c«(*) = MTriB^it) AU W (t)}] . (7) 

The procedure is repeated until the desired convergence has been reached. The hard numer- 
ical task is the propagation of the operators U and B with the time dependent Hamiltonian 



for which a second order Newton polynomial integrator [14 was used. 



A direct use of Eq. (0) in the algorithm leads to saturation. This is because the constraint 
related with the field energy becomes more important than the original objective. A remedy 
is to interpret the right hand side of Eq.(^), denoted by Ae^ k \ as a correction of the field in the 
previous interaction |T5|] . The field after the k iteration is then given by e^(i) = e^ _1 ^(t) + 
Ae (fc) (t). The OCT procedure was first applied to a 1 qubit operation the Hadamard rotation 
(Eq. |TD|). In a similar fashion, 2-qubit operation such as the "controlled not" were obtained. 
The optimal fields generating unitary operations with eight (3-qubits) were also attempted. 

It was found that the number of iterations required to converge the results increased 
approximately exponentially (or factorially) with the number of states in the problem. A 
possible reason for this scaling is that not only a specific state to state transition has to be 
forced but this has to be carried out without disturbing the other state to state transitions 
in U. As a result the scaling becomes 0(n\) which is consistent with numerical experience. 

In a molecular environment, obtaining the optimal field to carry out an algorithm is more 
involved. The register levels which are used to write the input and output are only part of a 
much larger manifold of molecular energy levels. Considering the advances in pulse shaping 
techniques in the visible region of the spectrum the transitions of choice are electronic. 
For such a molecular construction, imposing a unitary transformation on the total set of 
levels which are addressed by the field is too restrictive. Relying on the experience that the 
convergence is close to factorial, an extremely large number of iterations would be required 
to converge. 

The strategy is therefore to restrict the target objective to only the states used directly as 
registers keeping the other states in the system as passive observers. The reduced objective 
is obtained by changing the previous expressions Tr{AB} to ^(z|#AB|i) r , where is 
a basis of the subspace of registers. The subindex R is used to denote the operators in that 
subspace, for example Or denotes the target unitary transformation. The substitution of 
the restricted condition instead of the trace in Eq. (|j) keeps the maximum condition when 
the unitary transformation U#(T) is equal to Or. The maximum value becomes equal to 
the dimension of the subspace Nr. In the condition for B at time T, a particular dependence 
must be specified for all the levels. For simplicity the identity in the passive subspace has 
been chosen. 

The approach is illustrated using two models. For the first model, the implementation 
of a unitary transformation in one qubit while minimizing the population transfer to other 
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levels in the molecule is studied. The model consists of a molecule with two electronic 
surfaces described by the Hamiltonian, 



H= " . (8) 




where H s and H e are the ground and excited surface Hamiltonians. The electronic surfaces 
are coupled by the transition dipole operator p,, controlled by the shaped field e(t). The 
two first levels of the ground electronic surface are chosen as the registers representing the 
qubit. The model includes 15 rovibroinc levels in the ground electronic state and 5 in the 
excited state described by the Hamiltonians, 

15 5 

H e = ^2 E gi\9i)(9i\ ; H-e = ^2E ej \e j }(e j \. (9) 

i=l j=l 

shown in Fig. [I]. Next, a transition dipole operator with equal coupling strength between 
levels was chosen, p, = fi Yli=i Ylj=i ( e j\ + H.c, where H.c. denotes the Hermitian 
conjugate. 

In Fig. ^ the optimized field is shown when the target unitary transformation is a 
Hadamard rotation given by, 



(10) 



restricted to the 2x2 qubit subspace. The shape function was chosen as s(t) = sin 2 (27ri/T) 
and the initial guess was e^°\t) = s(t) cos(fit), where f2 is the frequency of the 00-line of the 
electronic transition. Fig. [I] shows that the dominant frequencies are the ones related to the 
known vibronic molecular transitions. In the time domain the field is split into a symmetric 
sequence of sub pulses. A phase relation correlating the dominant frequencies is observed in 
a Wigner plot (not shown). These phase relations guarantee that no population is lost to 
the excited levels. 

In the second model the two electronic surfaces, Eq. @, have two orthogonal vibrational 
modes, denoted as a and (3. The first two levels of each vibrational mode in the ground 
electronic surface are chosen as the the physical implementation of the two qubits. The 
goal is to induce an operation that involves entanglement between them. In this model 
each vibrational mode of the ground surface is coupled by the field to the corresponding 
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FIG. 1: (a) Schematic representation for executing the Hadamard rotation. Units are chosen such 
H = 1. The solid lines represent the two levels associated with the qubit. The levels in the ground 
surface are equally spaced with uj g = 1 (for Li2 in the A state uj g = 255 cm -1 ) and in the excited 
surface with u) e = 0.9 (for Li2 in the E state to e = 241 cm" 1 ). The other model parameters are 
Q = 15 (for h'i2 the A to E transitions f2 = 12400 cm -1 ) and = 0.1. The arrows are two of 
the possible transitions induced by the driving field, (b) Optimized field that induces a Hadamard 
rotation in the qubit at T = 70 (for L12 T ~ 1.5 psec). (c) Fourier transform of the field . 

mode in the excited electronic surface. The modes in the excited surface are coupled by 
a static term modeling Duschinsky rotation. This last term which is not controlled by the 
field can generate entanglement between the qubits. The model consists of two levels for 
each vibrational normal mode denoted by a and (3. The electronic surface Hamiltonians are 
H 9 = H ga <g> 1/3 + t a <g> H 9/3 and H e = H ea ® tp + l a <S> H e/3 + V Q/3 with, 

H 9i , = (E gUQ \g ) u (g \ v + E gui \g 1 ) u (g 1 \ u ) , 
H eu = (E euo \e ) u (e \ u + E evi \e\) v (e 1 \ v ) , 

V a/3 = 5 afS (|e ) Q ® |ei)^(ei| a ® (e |/j + |ei) a ® |e ) /3 (eoU ® (eil/3) , (11) 
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with V a /3 the Duschinsky term that couples the vibrational modes in the excited surface. 
A schematic representation of the levels is given in Fig. [|. The transition dipole operator 
is chosen as p, = il a <g> lp + l a <8> /Sp, with p, u = fi 0v {\e ) v + \ei) u )((g \„ + (gi\ u ) + H.c. 
Operators in the combined Hilbert space of this system are represented by 16 x 16 matrices. 




FIG. 2: (a) Schematic representation of the execution of a 2-qubit Fourier transform. The solid 
lines represent the levels associated with the qubits. The frequencies are E gao = 0, E gai = 1, 
Eea = 15, E eai = 15.8, E g p = 0, E g p 1 = 0.9, E e p = 14.5, and E e p x = 15.2. The arrows are 
examples of transitions induced by the driving field in each mode. The other model parameters 
are no a = 0.1, flop = 0.08 and 5 a p = 0.21. (b) Fourier transform of the optimal field that induces 

- FT 

O r at T = 320. 

The target unitary transformation is a two qubits quantum Fourier transform 0, 
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where the operator is represented in the basis \gi) a ® \9j)p- Notice that operators in the 
subspace of interest are represented by 4 x 4 matrices. The spectra of the field resulting 
from the optimization is shown in Fig. 0. The frequencies are confined in the region of 
the vibronic transitions. The structure of the spectra is much more complex than in the 
previous model. An obvious reason is the larger dimension of the subspace of interest. In 
addition the complicated mechanism required to create entanglement is indirect. During 
the pulse significant population is transfered to the excited electronic states. The resulting 
more complex field is reflected by also in the Wigner distribution, where the field shows a 
very high degree of correlation between time and frequency. 

In conclusion, Optimal Control Theory has been generalized to obtain the driving field 
that generates any target unitary transformation. In principle the scheme allows the im- 
plementation of any quantum gate using a single step. Compared to the implementation of 
the gate using a sequence of known simple pulses, the loss by the possible complexity of the 
shaped pulse can be more than compensated by the faster implementation. Convergence 
is efficient for a small number of levels. Besides, with the right choice of the guess field, 
the method can obtain the driving field for more complex systems where a large number 
of levels are involved and the pulses for simple operations are unknown, as in molecular 
systems. The advantage of such systems is the short time in which these algorithms can 
be executed. Using Li 2 as an example, the Hadamard rotation can be executed in ~ 1.5 
psec. The 2-qubit Fourier transform could be executed in a molecule for example OCS in 
approximately ~ 6 psec. These fast timescales give hope that the quantum computation 
can be carried out before decoherence processes take place. 
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